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SOME CONSIDERATIONS OF THE INFLUENCE OF BODY 
CROSS-SECTIONAL SHAPE ON THE LIFTING EFFICIENCY OF 
WING-BODY COMBINATIONS AT SUPERSONIC SPEEDS 
By E. B. Klunker and Keith C. Hander 

SIMMARY 


Linearized theory is used to estimate the order of magnitude of the 
beneficial effects that can be obtained by shaping the body of a wing- 
body configuration to generate lift indirectly at supersonic speeds. The 
analysis applies to a class of wing-body configurations which is believed 
to exhibit the essential features of this type of interference; the body 
is cylindrical upstream of a wing which has supersonic edges and the parts 
of the body above and below the wing are each semicircular and of different 
radii. Calculations are presented for a Mach number of |/~2 for a partic- 
ular family of bodies and a sonic-edge delta wing. 


INTRODUCTION 


At supersonic speeds, lift can be generated by locating a wing in 
the interference field of another component so that expansion waves occur 
over the upper surface of the wing and compression waves over the lower 
surface. Recent investigations have shown that in some cases efficient' 
lifting systems can be obtained by utilizing this interference. For 
example, in reference 1, Rossow considered the interference effects 
between a wing and planar elements which simulate a body or represent 
vanes at the wing tipB. Several combinations were determined which, 
according to linear theory, generate lift more efficiently t han the wing 
alone. However, as pointed out by Rossow, the friction drag of the fins 
may completely' offset the gain achieved. 

Aside from a wing-cone system analyzed by Rossow and a very limited 
amount of experimental work concerned with some effects of body cross- 
sectional shape, the studies dealing with the generation of indirect lift 
have been concerned primarily with the pressure fields and fordes devel- 
oped by interfering planar elements. Although the results of such studies 
are suggestive of the interference effects between wing and bodies, they 
do not apply directly to wing-body combinations wherein the body encloses 
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a certain volume and supports part of the lift. The purpose of the pres- 
ent investigation is to obtain some idea of the magni tude of the beneficial 
effects that can be obtained by shaping the body of a'ving-body configura- 
tion to generate lift indirectly. Accordingly, a configuration was selec- 
ted for study which is believed to exhibit the essential features of this 
type of interference and also yields to a relatively simple analysis. 

This configuration is composed of a wing with supersonic edges and a body 
which is cylindrical upstream of the wing, the parts of the body above 
and below the wing each being semicircular and of different radii. Cal- 
culations are presented for a Mach number of f2 for a particular family 
of body shapes and a sonic-edge delta wing.. 


SYMBOLS 


A(x) 

dA(x) 


b 

B 

c 

cl 


area distribution in upper half- space 



nondimensional wing span 
constant 

nondimensional wing chord 

lift coefficient based on total wing area 


Cj) drag'-coefficient based on total wing area 


D drag 

f(x) body source strength 

gpn(x) function defined in equation (6) 

h n (x) multipole strength 

L lift 

M free -stream Mach number 

p pressure 

q free-stream dynamic pressure 
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r nondimensional radius in. cylindrical coordinates 

Eo radius of cylinder 

r^ nondimensional body radius 

U free-stream velocity 

Wpn tabulated function 

x,y, z non dime nsional Cartesian coordinate system 

a wing angle of attack 

P = ju 2 - 1 

<t velocity potential 

0 disturbance velocity potential 

8 azimuth angle in cylindrical coordinates 

Subscripts: 

b refers to body potential 

1 refers to interference potential 

w refers to wing potential 

00 refers to free-stream conditions 

Superscripts : 
b refers to body 

v refers to wing 

ANALYSIS 

At supersonic speeds a lifting force on a wing-body configuration 
can be developed by contracting the body above the wing and expanding 
it below the wing. A configuration of this type which has the merit of 
simplicity of analysis consists of a flat-plate wing with supersonic 
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edges and a body which (l) is cylindrical upstream of - the wing and is 
at zero incidence, (2) is composed of semicircular parts of different 
radii above and below the wing, and (3) does not extend downstream of 
its intersection with the wing trailing edge. The simplicity of the 
analysis follows in part from the fact that the flow fields in the upper 
and lower half-spaces can be treated independently. The flow in each 
half-space can be represented by source distributions for the wing alone 
and body alone, together with an interference potential. With the body 
source distribution for 'one half-space the same as that for the other 
except for a change in sign, the body area at any station is equal to 
the area of the cylinder, and the pressure coefficients at correspo nding 
points above and below the wing are numerically equal but opposite in 
sign. Thus, the lift or drag for the part of the configuration in either 
half-space is one-half the total. 

The type of configuration con- 
sidered herein is shown in sketch a. 

The unit of length is the cylinder 
radius Eq; the quantities x, y, z, 

r = ^y 2 + z 2 } the wing chord c, and 

the span b are nondimens ional ; and 
the wing is at an angle of attack a. 

Since the analysis is similar to that 
for other wing-body problems treated 
extensively in the literature, only 
a brief outline of the formulation is 
presented here and the details of the 
calculations are relegated to the 
appendix. The problem is formulated 

for the upper half- space. Sketch a 



The velocity potential $ is taken as 


® = URo(x + 0 W + 0 b + 0i) 


( 1 ) 


where the disturbance potentials satisfy the linearized differential 
equation for the velocity potential, 0^. and 0 b are the potentials 

for the wing alone and body alone, respectively, and 0^ is an inter- 
ference potential required to satisfy the condition of zero normal veloc- 
ity at the surface of the wing-body combination. The requirement of zero 

normal velocity on the wing is ~ = -a With the definitions of the 

oz ax 
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disturbance potentials together with the fact that 

zero for z = 0 for a circular body, this boundary 
mated by 


btfbjbz is identically 
condition is approxi- 


3z 


-a 



(z = 0 ) ( 2 ) 


3<f 


dr> 


3<I> 


The requirement of zero normal velocity on the body is — 

ox dx ' or 

at r = r^,, where is the nondimens ional body radius. With the 

restriction that the body does not depart substantially from the cylinder, 
this boundary condition is approximated by 


a^ + aA =0 (r . x) (5) 

Or dx Or Or 


The body potential is represented by a distribution of sources along 
the body axis, and in the linear approximation the pressure p^ due to 

the body is 


Pfc> “ P<» 

” i " 


1 f x_pr f(i) cli 

0 \j(x - Z ) 2 - fa 2 


W 


where the prime denotes the derivative with respect to the indicated 
argument, [3 = ^M 2 - 1, and M is the stream Mach number. The source 
strength f(|) is related to the body geometry by the first of equa- 
tions (3). An equivalent statement of this boundary condition for quasi- 
cylindrical flows relates the body area in the upper half- space A(x) to 
the source strength by 


A(x) - ~~= AA(x) = J ^ f’(l)^(x - t ) 2 - p 


cL| 


2 


(5) 
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where AA.(x) is the change in area in the upper half-space. Since the 
body source strength for the flows in the upper and lower half-spaces 
differ only in sign, AA.(x) for each half-space is the same except for 
a change in sign, and the total body area at any station is equal to the 

O 

area of the cylinder jtRq . 

The interference potential is required tQ cancel the flow through 
the body due to the wing potential and is represented by a distribution 
of multipoles along the body axis as 



f 


cos 


i^0t sin 


n8 


f 0 


x-pr IfeU) 



where h n ( |) is the multipole strength and 0 is the azimuth angle 
defined in sketch a. The boundary condition on as given by equa- 


tion (3) is satisfied in each half-space by expanding 


5 0w 

Sr 


in a Fourier 


series and equating the coefficients of each harmonic to the corresponding 

coefficient of . Because of the symmetry properties of the wing, only 

Sr 

the cos 2n0 terms appear, and from equation (3) 


- - — = - — = a 2_ S 2 n( x ) cos 2110 
Sr Sr n=:0 


(r - 1) (6) 


where the Fourier coefficients are given by 




d0 


r=l 




cos 2n0 d0 


r=l 


t* 
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The interference potential with <j>^ « cos 2n8 satisfies the boundary- 
condition as given by equation (2) for (j>± identically. The multipole 
strengths can be deter mine d from equation (6) by replacing by the 

expression for a distribution of multipoles. However, this step is not 
necessary since only the pressure is required and Nielsen (ref. 2) has 
shown that the interference pressure for a quasi- cylindrical body can be 
related directly to the functions gpn(x) and a tabulated function 

Wpn(x,r) • This expression for the interference pressure p.j_ is 



— 2_ cos 2n0 -i; g^x-pr+p) - 

P r£0 ^ 


1 

P 


^-p(r-l) 


P 



r+l,rj 


d| 


( 7 ) 


The pressure p w due to the wing is determined directly from the stream- 
wise disturbance velocity of the wing alone. 

The lift and drag are determined by integrating the surface pres- 
sures. The total lift L is 


L 

2Rq 2 


'II ( P t + P v + p l) “ • 

Wing 

panels 




L b 

2EQ 2 


( 8 ) 


where the superscript on the pressure refers to the surface where the 
pressure acts, w denoting the wing and b the body. The superscript 
(or b) on a force refers to the force obtained by integrating over the 
wing (or body) . Similarly the total drag D is 
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D 

2Eo 2 


a L w 1 

=—r + - 
2% 2 * 



ZA'(x) dx 



P§ + + P 


d0 


D y + _D^_ 
2Eq 2 2Eq2 


(9) 


■where the first term is the drag of the wing panels and the second term 
is the body drag. 


DISCUSSION 


The method presented in the preceding section affords a means for 
evaluating the lift and drag characteristics of a class of -wing-body 
configurations. For a given wing plan form, the forces can be determined 
most easily by specifying the body source distribution rather than the 
body geometry. Calculations have been made for a Mach number of ^2 for 
a delta wing with sonic leading edges and a linear body source distribu- 
tion. This configuration was selected primarily to simplify the calcula- 
tions. However, it is believed that the results are ' generally indicative 
of the effect of a body cross-sectional shape on the lifting efficiency of 
a wing-body configuration at supersonic speeds. 


The body area distribution is determined from equation ( 5 ) for the 
body source distribution given by f(x) = -Bx, where B Is a positive 
constant. The body shape Is shown in 
sketch b for the rather large value of 
B = 0-55* The contributions to the 
lift and drag from the disturbances, 
arising from the wing and body can be 
integrated directly and are given in 
the appendix, whereas the contribution 
to the forces from the interference 
pressure must be determined numerically 
The interference pressure depends only 
upon the geometric properties of the 
wing through the Fourier coefficients 
g 2n ( x ) . These functions have been 

evaluated for n = 0 and 1 in refer- 



ence 2 and only the first two terms of the interference pressure have 
been used in the calculations. Although this limitation affects the 
accuracy of the calculations, Nielsen indicates that the major effect 
is contained in the first two terms. 
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The lift and drag coefficients, based on wing area, can be written 
as the sums 


C L = °L,w + c L,b + °L,i + °L,w + c L,b + c L,i 

c D = -(ct,: + C L J + C L j) + + Oo f l + C D ^ 


where the subscripts denote the origin of the pressure disturbance and 
the superscripts denote the surface upon which the pressures act. The 
lift coefficients with the subscripts w and i are proportional 

to a and those with subscript b are proportional to B, while ^ 

and 8X6 Proportional to Ba and Cp ^ is proportional to B^. 

Consequently, the ratio Cp Jcj^, which is indicative of the lifting 

efficiency of the configuration, is a function only of B/a and c. 

The values of B/a which correspond to the minimum values of 

are denoted by (B/a)^^ and are shown in figure 1 aB a function of c. 
In figure 2, values of Cq / cp, 2 are presented for the configuration with 
(B/a) 0 p^. together with those for the wing alone, the wing-cylinder 

(B = 0), and the indented cylinder with wing at zero angle of attack 
(a = 0). The curve for the wing with the optimum Indented cylinder is 
the locus of points for body shapes which give rise to the least values 

of thus the body shape varies with angle of attack and c in 

accordance with figure 1. The curve for the wing alone, which exhibits 
the highest lifting efficiency, is presented primarily for the purpose 
of orientation - the most meaningful comparisons are those between the 
configurations which have a body. 

Comparison of the curves for the wing-cylinder (B = 0) and the 
indented cylinder with wing at a = 0 shows that for large wings (that 
is, large values of c) lift is generated more efficiently by the wing- 
cylinder, whereas for small wings lift is generated more efficiently by 
shaping the body. Limiting calculations for large values of c show 
that the lifting efficiency of the wing-cylinder approaches that of the 
wing alone, while for the indented cylinder Cp j approaches infinity 

as c approaches infinity. The lifting efficiency that can be obtained 
by generating lift both directly and indirectly represents a sizable 
gain over the efficiency obtained with either method individually. The 
decrease of (B/a) 0 p t with increasing c reflects the fact that, for 
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the optimum condition, proportionately more of the lift is developed by 
the system which itself is most efficient. 

The quantities Cp^x, and Cp/a 2 for the wing, the wing- cylinder, 
and the configuration with (B/a) Q p^. are presented in figure 3 as func- 
tions of c. The value of Cp^a for the configuration with (B/a) 0 p^ 

is substantially greater than for the wing-cylinder and is nearly equal 
to that of the wing alone for c equal to 5* The quantity Cp/a 2 for 

the configuration with (B/a) 0 p^ is approximately equal to that for the 

wing alone for c = 2, and becomes approximately 10 percent higher at 
the larger values of c. The value of Cp^a 2 for the wing- cylinder is 

substantially lower since the cylindrical body does not contribute to 
the drag. The differences between the curves of figure 3 for the con- 
figuration with (B/a) opt and the wing-cylinder represents the contri- 
bution to the forces due to shaping the body. 

The contributions to the lift and drag coefficients are given in 
figures 4 and 5 for the optimum configuration. The contributions to the 
lift and drag from the interference pressure produces a down load and 
thrust. It is of interest to note that, for the optimum condition, the 

v b v 

lift coefficient Is closely equal to Cp w + b since "the terms 1 

and Cp^p nearly cancel Cp^p and Cp^£, respectively. Similarly, for 

V V 

the optimum condition, the terms and C-p^ p nearly cancel while 

v j_ 

Cp^ w and Cp^p together contribute a small value which Is nearly inde- 

b b 

pendent of c. The terms Cp^ b and decrease in value with 

increasing c, reflecting the fact that for the optimum condition B/a 
decreases as the ratio of wing span to cylinder radius increases. 

The effectiveness of using indirect lift to increase the lifting 
efficiency of wing-body comb ina tions is dependent to a large extent on 
the body configuration. Although only the portion of the body in the 
neighborhood of the wing contributes to the indirect lift, the shape 
ahead of the wing affects the body drag and consequently the optimum 
ratio of direct to indirect lift. Nevertheless, the example presented 
serves to illustrate the effect of body cross-sectional shape on the 
lifting efficiency. The family of body shapes used in the example was 
selected for simplicity of analysis. It seems clear, however, that other 
area distributions would be superior for developing indirect lift. For 
example, the body shape in the neighborhood of the wing (for each half- 
space) corresponding to a square-root type of source distribution, a type 
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which arises in various minimum-drag problems for bodies, would probably 
lead to a lifting system more efficient than that of the example since 
it would be favorable to both the body drag and the wing lift. 

The effectiveness of developing lift by utilizing the pressure field 
from the body is dependent upon the wing plan form and upon the stream 
Mach number. The increase in lifting efficiency obtained by utilizing 
indirect lift might be expected to be somewhat higher than that of the 
example for wing plan forms, such as swept wings, which could make use of 
more of the pressure field arising from the body to develop lift. More- 
over, it would appear that body shaping is most effective for the con- 
dition where the wing leading edges are sonic, since the pressure relief 
for wings with subsonic edges and the smaller wing area affected by the 
pressure field from the body for wings with supersonic edges would result 
in less lift from the interfering body elements. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , June 18 , 195 &. 
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APPENDIX 

EVALUATION OF TEE PEESStIRES AND FORCES FOR THE EXAMPLE 


The expressions for the pressures and forces for a delta wing on a 
quasi- cylindrical body are presented herein. The body shape for each 
half-space corresponds to that given by a linear source distribution. 

The calculations are for the sonic edge condition and a Mach number of \[ 2 . 

From equations (8) and ( 9 ) the lift and drag coefficients for this 
configuration are 


01 ■ " li* fi * 4° * p; * i? ) “ ■ j s° ^ r ( p ‘ + * 5 ») “* 8 “ 

/ y ° (»; . p; * ^ ^ 4° *•<«> « 4’ (pg * pS <• p|) « 


(Al) 


For a linear source strength f(|) = -B|, where B is a positive 
constant, the body shape in the upper half-space is determined from equa- 
tion (5) as 


ZA(x) 



(A2) 


The pressure due to the body is determined from equation (4) as 



~ cosh"’*' 

Tt 


X 

r 


(A3) 


and and p^ are given by equation (A3) for r = 1 

respectively. 


and r = y. 
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The line pressure sources of reference 3 can be used to determine 
the pressure due to the wing. For the sonic edge condition and = 1 , 
Py. becomes 


Pw ~ Poo 
q " 



W 


From equation (A^), p^ is obtained for r = 1 and p^ is obtained 
for 0 = 0 and r = y. 

The interference pressure must be evaluated numerically since it 
involves the tabulated functions g 2 n.( x ) an<i Wpy,. These functions 

have been evaluated for a flat-plate wing in reference 1 for n = 0 
and n = 1 . The functions g 2 n( x ) used herein correspond to the func- 
tions -f pna (x) given in figure 2 of the reference; the function gg(x) 
is a constant equal to 2 /jt. The interference pressure on the wing p^ 

is determined from equation (7) by retaining only the first two terms of 
the infinite series for p.|_, and setting z equal to 0 and 0 equal 
to 0: 



_1 

>_ 

n=0 


l , nx-y+l 

4 = ggnU-y+l) - J S2n^) w 2n( x -6-y+l^y) < 3-1 




(A3)- 


The pressure on the body 


„b 

p i 


is given for r = 1 as 


b 

Pi " Pa 


= 2a 2_ cos 2n0 
n=0 


g2n( x ) - ggaCO WgnCx-l,!) d| 


(AS) 



Ik 
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When the expressions for the body shape and pressures are substituted 
into equations ( Al ) , the various contributions to the lift and drag of 
the wing are 



and the contributions to the lift and drag of the body are 



2 

? 

qac 




sin 0 d9 = 



cosh' 


c - 






p£ sin 9 d9 



coth“-*- 



cosh 


-1 
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where 

Kl) = 

V. 

2 

a 

2 

qac^ 

r b 
°D,b 

1 

a 2 

cjn:c 2 a, 2 





0 b 
l D,w 

i 

a 2 

qrtcSa 2 ■ 


— -U 


itc 2 a V 

°D,i 

1 

a 2 

qrtc 2 a, 2 1 


3*5' 3*5*5* 3*7*7* 3-5-9J 


fax f 

Jo 


__ k-/ 

Pi 


sin 0 d.9 


cosh' 


' 1 c^2c \l< 


'c* - 1 - cosh - " 1 ' c) + 1 - c 


- — r AA’(x) dx ( p b d0 

: 2 a 2Jl Jo 1 


ire 


ii_ B 

2 a| 


i^c 2 - 1 - cosh -1 c^ - ^rn dx J' ¥ 0 (x-i,l) d| 


The coefficients for 2b 1 ., and 25ii 

cty a 2 y a a 2 


have been integrated 


numerically by vising the first two terms of the series representation. 
The various contributions to the lift and drag are presented in figures 
and 5 as a function of c i^^JiaBt:ws3aeaKffbf B/a, which mi nimize Cj 


ro4=- 
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